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^ ' Abstract 

We revisit proton decay via the color-triplet Higgs multiplets in the minimal 
supersymmetric grand unified model with heavy sfermions. Although the model 
has been believed to be excluded due to the too short lifetime of proton, we have 
found that it is possible to evade the experimental constraints on the proton decay 
rate if the supersymmetric particles have masses much heavier than the electroweak 
—^ ' scale. With such heavy sfermions, the 126 GeV Higgs boson is naturally explained, 

while they do not spoil the gauge coupling unification and the existence of dark 
matter candidates. Since the resultant proton lifetime lies in the regions which may 
be reached in the future experiments, proton decay searches may give us a chance 
to verify the scenario as well as the supersymmetric grand unified models. 
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1 Introduction 

The discovery of the Higgs boson [T||2] has opened the way for physics beyond the Standard 
Model (SM). It is certainly a striking hint for understanding the high-energy physics, as 
elementary scalar particles may play an important role in realizing our complicated world 
with apparent broken symmetries. The theories with supersymmetry (SUSY) naturally 
include such scalar particles; the Higgs boson might be a superpartner of chiral fermions, 
called higgsinos. Besides, there exists a set of scalar particles for each SM fermion, as 
well as adjoint fermions for the SM gauge bosons. Then, astonishingly enough, we find 
that with these extra particles the gauge coupling constants of the SM are to be unified 
at a certain high-energy scale with great accuracy [3_H7]. This observation motivates us 
to study the supersymmetric grand unified theories (SUSY GUTs) [SJE]. 

The SUSY GUTs predict an exciting phenomenon: proton decay. It is induced by 
the exchanges of the color-triplet Higgs multiplets and the X-bosons, and their effects 
are described in terms of the dimension-five and -six effective operators, respectively. In 
the minimal SUSY SU(5) GUT [HIE], which is a simple supersymmetric extension of the 
original SU(5) GUT [10], the former process yields the dominant decay modes, such as 
p — > K + v. The lifetime of the channel is estimated as r(p — > K + v) < 10 30 yrs [TTjITS] . with 
the SUSY particles, in particular those of the third generation, assumed to have masses of 
around the electroweak scale. On the other hand, the Super-Kamiokande experiment gives 
stringent limits on the channels: r(p -> K + v) > 4.0 x 10 33 yrs [13J. This contradiction 
makes it widely believed that the minimal SUSY SU(5) GUT has been already excluded 
and, therefore, needs some extensions in order to suppress the dimension-five proton decay. 

As is often the case with SUSY models, the SUSY GUTs are usually discussed within 
the context of the low-scale supersymmetry. Recently, experiments at the Large Hadron 
Collider (LHC) provide limits on the SUSY models. The ATLAS and CMS Collaborations 
have been searching for the SUSY particles and imposed severe constraints on their masses, 
especially those of squarks and gluino [T4T - T16] . The mass bounds have began to exceed 
1 TeV and, thus, the low-energy SUSY models are confronted with difficulties. Moreover, 
the observed mass of the Higgs boson around 126 GeV [U|2] might also indicate the SUSY 
scale is considerably higher than the electroweak scale; in the minimal SUSY Standard 
Model (MSSM), the mass of the lightest Higgs boson is below the Z-boson mass at tree 
level, so sufficient mass difference between stops and top quark is required in order to 
raise the Higgs boson mass through the radiative corrections [TTH2T| . 

In fact, the SUSY models with heavy SUSY particles have a lot of attractive features 
[22TI28] . First of all, since the flavor changing neutral current (FCNC) processes and/or 
the electric dipole moments induced by SUSY particles are suppressed by their masses, the 
SUSY flavor and CP problems [29J are relaxed when the masses are considerably heavy. 
As mentioned to above, the heavy sfermions yield sufficient radiative corrections to lift the 
Higgs mass up to 126 GeV [30H33] . They do not spoil the gauge coupling unification since 
the sfermions form complete SU(5) multiplets. Actually, it turns out that the unification is 
improved in the sense that the required threshold corrections at the GUT scale tend to be 
reduced [34J . As for the cosmology, the gravitino problem is avoided because of the high- 



scale SUSY breaking, and the thermal leptogenesis scenario well works with high reheating 
temperature [32]. Further, this high-scale SUSY scenario naturally accommodates dark 
matter (DM) candidates, which might be detected in future dark matter experiments 
directly [36H39] and indirectly [32JH0JI1T] . Thus, with the recent LHC results considered, 
the high-scale SUSY scenario is even promising from a phenomenological point of view. 

Interestingly, this scenario also provides an alternative solution to the problem re- 
garding the dimension-five proton decay in the minimal SUSY GUT. The dimension-five 
operators generated via the color-triplet Higgs exchange contain squarks and/or sleptons 
in their external lines. These fields are to be integrated out below the SUSY scale through 
the wino or higgsino exchanging processes, and then the four- Fermi operators, suppressed 
by the sfermion masses, are induced. Hence, their effects are expected to be extremely 
reduced when the SUSY scale is much higher than the electroweak scale. 

In this paper, we study such possibilities within the context of the high-scale SUSY 
scenario. We will find that the minimal SUSY SU(5) GUT actually evades the constraints 
from the proton decay experiments with the SUSY braking scales which naturally explain 
the 126 GeV Higgs boson and the existence of dark matter in the Universe. The resul- 
tant proton lifetime lies in the regions which may be reached in the future proton decay 
experiments. Therefore, although the high-scale SUSY scenario is hard to be probed in 
the collider experiments, the proton decay searches may give us a chance to verify the 
scenario as well as the existence of supersymmetry and the grand unification. 

This paper is organized as follows: In Sec. El a high-scale SUSY model which we discuss 
in this work and its phenomenology are briefly explained. In the next section, we give 
a set of formulae for evaluating the proton decay rate via the dimension-five operators. 
Then, we show the resultant proton lifetime in the model, and compare it with current 
experimental limits in Sec. HI Section [5] is devoted to conclusions and discussion. 

2 High-scale SUSY 

To begin with, we describe a high-scale SUSY model which we deal with in the following 
discussion. Assume that there exists a SUSY breaking hidden sector where the SUSY 
breaking is triggered by a chiral superfield Z which is not a gauge singlet. Then, with 
a generic form of Kahler potential, all the scalar bosons except the lightest Higgs boson 

acquire masses of 

F 7 

M s , (1) 

with Fz and M* the F-component vacuum expectation value (VEV) of the field Z and 
the mediation scale of the SUSY breaking, respectively. The gravitino mass is 7713/2 = 
Fz/v3M-p\, and it is the same order as the scalar masses when M* is around the Planck 
scale, M-p\. The gaugino masses, on the other hand, are not generated by the dimension- 
five operators like 

r d 2 9 ZTi[W a W a ] , (2) 



since the symmetry under which the superfield Z is charged prohibits such an operator. 
Here, W a denotes the gauge field-strength chiral superfield. Instead, they are induced by 
the anomaly mediation mechanism [4*2"||4"3] : 



M a = b ^Lm 3/2 , (3) 

where Mi, M 2 , and M 3 are the masses of bino, wino, and gluino, respectively, and b a are 
the one- loop beta- function coefficients of the gauge coupling constants g a (a — 1, 2, and 
3 for U(l)y, SU(2)l, and SU(3)c, respectively). This expression tells us that the gaugino 
masses are suppressed by one-loop factors compared with the gravitino mass. Finally, the 
higgsino mass, \xu, is somewhat model-dependent; it might lie around the same order of 
gaugino masses when some additional symmetries exist, or be as large as gravitino masses 
if it is generated by the Kahler potential. Thus we regard it as a free parameter in the 
following discussion. 

In this scenario, the lightest SUSY particle (LSP) is either wino, which turns out to 
be the lightest gaugino, or the lighter higgsino. Then, it is found that in any case the 
LSP may explain the dark matter in the Universe. Indeed, the thermal relic abundance 
of wino and higgsino DM with a mass of 2.7-3.0 TeV [44] and 1 TeV [45], respectively, 
accounts for the observed density of DM. With relatively small masses, the non-thermal 
production of them also might be consistent with the observation [4*61147] . The wino or 
higgsino DM in this model implies sfermion masses lie around 10 2 -10 3 TeV. 

From now on, we assume the sfermions are nearly degenerate in mass, and their masses 
are collectively denoted by Mg. The mass is supposed to be Mg — 10 2 -10 3 TeV, and either 
wino or higgsino is assumed to be the LSP. Models with such a mass spectrum and their 
phenomenology have been enthusiastically studied in the previous literature [451455] . 

As mentioned to in the Introduction, the mass spectrum does not spoil the gauge 
coupling unification [23j[2l], since the sfermions are embedded in the complete multiplets 
of SU(5). In fact, the unification may be improved [31]; a renormalization group analysis 
reveals that all of the GUT scale particles, especially the color-triplet Higgs multiplets, 
possibly lie around the GUT scale ~ 10 16 GeV, contrary to the case of the low-energy 
SUSY. It indicates that the threshold corrections to the gauge coupling constants at the 
GUT scale are reduced in the high-scale SUSY scenario. Thus, the SUSY GUTs are still 
well-motivated. 

3 Proton Decay in the Minimal SUSY SU(5) GUT 

In this section, we review the proton decay in the minimal SUSY SU(5) GUT. In this 
model, the MSSM matter fields are embedded in a 5 © 10 representation. The SU(2)^ 
singlet down-type quarks D and doublet leptons L are incorporated into the 5 fields, $, 
while the SU(2) L singlet up-type quarks, U, doublet quarks, Q, and singlet leptons, E, 
are formed into the 10 representations, \1/. Here all the superfields are expressed in terms 



of the left-handed chiral superfields. The explicit form of the multiplets is 
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Here, a = 1,2,3 denotes the color index. The MSSM Higgs superfields, on the other 
hand, are embedded into a pair of 5 and 5 superfields accompanied with the new Higgs 
superfields Hq and Hca called the color-triplet Higgs multiplets: 
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where the last two components are corresponding to the MSSM Higgs superfields, 

H„ I 7,o ) • H d = (S) . (7) 



Exchanges of the color-triplet Higgs multiplets induce the baryon-number violating 
interactions. They are coupled with the ordinary matter fields by the Yukawa coupling 
terms in the superpotentiao 



W- 



Yukawa 



-h ij e abcde ^f V? H e - V2f ij ^% a H b 



(8) 



where i,j = 1, 2, 3 and a,b,c, . . . = 1-5 represent the generations and the SU(5) indices, 
respectively. The Yukawa couplings h l i and / y in Eq. (jHJ) have redundant degrees of 
freedom, most of which are eliminated by the field re-definition of \1/ and $ [56j. We 
parametrize the couplings according to Ref. [57] as 
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(9) 



1 In this paper we evaluate the proton decay rate with the SU(5) symmetric Yukawa couplings which 
are evaluated with up-type and down-type quark masses and the Cabibbo-Kobayashi-Maskawa matrix, 
though the ratios of charged lepton and down-type quark masses are not necessarily consistent with them. 
We have checked that, even if the lepton masses are used for fd i7 our consequence presented below is not 
changed significantly. 
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Figure 1: Supergraphs which illustrate color-triplet Higgs exchanging processes where 
dimension-five effective operators for proton decay are induced. Bullets indicate color- 
triplet Higgs mass term. 



with Vij the Cabibbo-Kobayashi-Maskawa (CKM) matrix. The phase factors ipi are sub- 
ject to a constraint 

<fi + ¥2 + y?3 = , (10) 

and thus two of them are independent parameters. The Yukawa coupling terms in Eq. (JHJ) 
are written with the component fields as follows: 

^Yukawa = Jm^^rJJ ioQT H^ ~ V*jfd j £rsQ'i' r Dj a Hj — fd j e rs V*jEiL'jH i J 



+ f^UiaEtHg - V£f di e c WU ia D j f i Hci • 



-f u ^e aM e rs QTQTHl + Vtf dl e rB QrL' i H a a 



l] J Uj 



(11) 



Here, r, s are the SU(2) L indices. In the following discussion, we also use the flavor basis 
for the matter fields. The relations between the flavor and gauge eigenstates are given as 
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(12) 



where primes represent the flavor eigenstates. 

The Yukawa interactions of color-triplet Higgs multiplets give rise to the dimension- 
five baryon-number violating operators. The processes in which the operators are induced 
are illustrated by the diagrams in Fig. [TJ After the color-triplet Higgs multiplets are 
decoupled, we obtain the effective superpotential as 
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Figure 2: One-loop diagrams which yield the baryon-number violating four-Fermi op- 
erators. Diagrams (a) and (b) are generated by charged wino and higgsino exchanging 
processes, respectively. Gray dots indicate dimension-five effective interactions, while 
black dots represent wino or higgsino mass terms. 



which yields the dimension-five effective operators, 



£ 5 = / d 2 9W 5 + h.c. 



(14) 



The effective operators contain sfermions in their external lines. Below the SUSY 
breaking scale, Ms, these sfermions turn into the SM fermions via the charged wino and 
higgsino exchanging processes shown in Fig. [2] In this figure, the gray and black dots 
indicate the dimension-five effective interactions and the mass terms for wino or higgsino, 
respectively. The first operator in Eq. (fT3j) contributes to the diagram (a), while the 
second one induces the diagram (b). Although the contribution of the diagram (b) is 
suppressed by the CKM matrix elements as it is generated in the flavor changing process, 
it is found to be sizable because of the large Yukawa couplings of the third generation 
fermions [H1EB]- The contributions of flavor-conserving neutral gauginos and higgsino 
exchange are in general suppressed by the Yukawa couplings of the first generation, thus 
negligible. Among them, the gluino contribution might be sizable because of the large 
coupling. It turns out, however, that the gluino contribution vanishes in the limit where 
squarks are degenerate in mass, and we consider such a case in the following calculation. 
When sfermion mass matrices have large flavor mixing, the contributions may also be 
significant, though we do not take into account such a situation for simplicity. 

After the electroweak symmetry breaking, the charged wino and higgsino are mixed 
with each other. In the following calculation, however, we neglect the effect since we 
mainly consider the case where M2, [in 3> f^w with mw the mass of H^-boson. When the 
masses of wino and higgsino are nearly degenerate, the mixing effects might be significant. 
It is straight-forward to modify the formulae obtained below in such a case. 



Evolving the four-Fermi operators from the SUSY breaking scale to the hadron scale 
(~ 1 GeV) according to the renormalization group equations (RGEs), we finally obtain 
the effective operators for proton decay as 
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where we use the two-component spinor notation for the SM fermion fields; all of the 
quarks are written in the flavor basis though primes are omitted for simplicity; Mg is the 
mass of sfermions with all the sfermions assumed to be degenerate in mass; m q are the 
masses of quarks defined in the DR scheme at the scale of /i = 2 GeV; u 2 , d 2 , u 3 , and d 3 
denote c, s, t, and b quarks, respectively; a 2 = g 2 /^ with g 2 the SU(2) L gauge coupling 
constant at the electroweak scale; tan/3 = (H^)/(H^); A£ and Ar in Eq. ( TTBT) represent 
the renormalization factors. These factors include the renormalization effects for both 
the couplings and the effective operators. The estimation of these factors is carried out 
in Appendix lAl 

The loop function F(M, Mg) is given as 



F(M, M 2 S ) = M 



1 M 2 , fM\ 



In 
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(16) 



where M is either the wino mass, M 2 , or the higgsino mass, [in- In the limit of M <^ Mg, 
the function leads to 

F(M,M|)^^, (M iS ^>M), (17) 

while in the limit of M — > Mg, it follows that 

F(M,MJ)^^-, (M^Mg). (18) 

Note that the function is proportional to M, when M ^ Mg. For this reason, the contri- 
bution of the diagrams in Fig. |2] is enhanced when the masses of the exchanged particles 
are large. In particular, in the case of fin ^ M 2 , the higgsino exchange contribution 
(the diagram (b) in Fig. [2]) dominates the wino exchange one. We also find from the 
behavior of the loop function that the transition amplitude is considerably suppressed 
when the sfermions have sufficiently large masses. Thus, we expect that the experimental 
constraints on the proton decay rate may be avoided in the high-scale SUSY scenario. 

The effective operators in Eq. ( TT5|) are written in terms of partons. In order to derive 
the decay amplitude for proton, we need to obtain the matrix elements of the quarks 



appearing in the operators between the proton state \p) and the kaon state \K + ). With 
the matrix elements, we finally obtain the partial decay widths, Y(jp — V K + P fi ) and F(p — >■ 
K + u T ), and the sum of them well approximates the total decay rate of the p — > K + v 
channel. Explicit formulae for the decay widths as well as the determination of the 
matrix elements are presented in Appendix. Q3] 

Before concluding this section, we briefly comment on the proton decay via the SU(5) 
gauge boson exchange. The SU(5) gauge bosons, called X-bosons, give rise to the 
dimension-six operators which contribute to proton decay by the p — > 7r°e + channel. In 
Ref . [34J , it is pointed out that the GUT scale in the high-scale SUSY tends to be slightly 
lower than that in the low-energy SUSY. Thus, the proton decay rate in this channel is ex- 
pected to be enhanced. It turns out, however, that the resultant lifetime is generally long 
enough [SH] to evade the current experimental bound, r(p — > 7r°e + ) > 1.29 x 10 34 yrs [50] . 
Thus, we ignore the contribution in the following calculation. 

4 Results 

Now we show numerical results of the proton decay lifetime in the high-scale SUSY sce- 
nario. We will see below that the resultant lifetime is well above the current experimental 
limits in a wide range of parameter region. 

First, we consider the case where the higgsino mass is of the order of the sfermion 
masses, Ms- In this case, the higgsino exchange contribution (the diagram (b) in Fig. [2]) 
dominates the wino exchange one, as mentioned to in the previous section. For this 
reason, the lifetime has little dependence on the additional phases, (fi in Eq. (Q, as well 
as the wino mass. Thus, it is possible to make a robust prediction for the proton decay 
lifetime. As the right-handed stop and stau run in the loop in the higgsino exchanging 
diagram, Ms should be regarded as their masses, which we assume to be degenerate for 
brevity. For Ms = [in-, the proton lifetime t p is approximately given asj 



^4x 10 35 x sin 4 2/3 ^ -±* — -£— yrs , (19) 



,35 „ ^4o« / 0-lVf M S V M Hc 

A R ) V!0 2 TeVy VlO 16 GeV 

and found to be well above the current experimental limits, r(p — > K + v) > 4.0 x 10 33 yrs 
[15] . with the SUSY scale being much higher than the electroweak scale. 

Let us investigate it in detail. To begin with, we consider the mass of the color-triplet 
Higgs multiplets, Mh c - As mentioned to in Sec. [2j through the RGE analysis discussed 
in Refs. [5T]l6T] with requiring the gauge coupling unification, one finds that Mh c may be 
around the GUT scale in the case of high-scale SUSY [51]. The prediction is, however, 
quite sensitive to the mass spectrum below the GUT scale, especially to the masses of 
higgsinos and gauginos. Thus, in the following discussion, we just fix M Hc to be around 
the GUT scale. One easily obtains proton lifetimes corresponding to other values of Mh c 
by using the power law given in Eq. f fl9|) . 



2 The renormalization factor reduces the proton decay rate for larger Mg when M$ = Hh = (10 2 - 
10 5 ) TeV, as described in Pig. [5] in Appendix. [A] 
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Figure 3: Lifetime of p — >■ f^ + z/ mode as functions of M5 = /i#. Wino mass is set to 
be 3 TeV and M Hc = 1.0 x 10 16 GeV. Solid lines correspond to tan/3 = 3, 5, 10, 30, and 
50 from left-top to right-bottom, respectively. Shaded region is excluded by the current 
experimental bound, r(p — > K + v) > 4.0 x 10 33 yrs |13j . 



In Fig. [31 we present the lifetime of the p — > K + v mode as functions of Ms = ^h- Here, 
the wino mass is set to be 3 TeV, while the result scarcely depends on the mass as long as 
M 2 < M s . The color-triplet Higgs mass is fixed to M Hc = 1.0 x 10 16 GeV. The solid lines 
are for tan/3 = 3, 5, 10, 30, and 50 from left-top to right-bottom, respectively. The shaded 
region is excluded by the current experimental bound, r(p — > K + v) > 4.0 x 10 33 yrs |13j . 
The figure illustrates the behavior presented in Eq. (TT9l . Moreover, it is found that 
the proton decay lifetime in the high-scale SUSY scenario may evade the experimental 
constraints, especially for small tan/3 and high SUSY breaking scales. We also show a 
similar plot for a relatively small value of Mh c in Fig. HI where the mass is taken to 
be 1.0 x 10 15 GeV. The wino mass is again set to be M 2 = 3 TeV, and the solid lines 
correspond to tan/3 = 3, 5, 10, and 30 from left-top to right-bottom, respectively. We see 
that the relation between the results presented in Figs. [3] and H] is well explained by the 
simple power law in Eq. f|T9|) . though the renormalization factors may also be changed 
with different values of Mn c . For this reason, we just fix Mh c = 1.0 x 10 16 GeV in the 
following analysis. One easily read other results with different values of Mh c by using 
the relation given in Eq. ( 1T9|) . 

Next, we consider the case where the higgsinos are lighter than the sfermions. In 
this case, the lifetime depends on the new phases appearing in Eq. ([9j). Here, we take 
the phases so that they yield the maximal amplitude for the proton decay rate, i.e., we 
require that each term in Eqs. (148 j) and (149 j) be constructive. This requirement together 
with the constraint ([TO]) uniquely determines all of the phases tpi. Since the choice of 
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Figure 4: Lifetime of p — > K + v mode as functions of Ms = fin- Wino mass is set to 
be 3 TeV and M Hc = 1.0 x 10 15 GeV. Solid lines correspond to tan/3 = 3,5,10, and 
30 from left-top to right-bottom, respectively. Shaded region is excluded by the current 
experimental bound, r(p — > K + v) > 4.0 x 10 33 yrs [T3] . 



phases gives the maximal proton decay rate, we are to obtain the most stringent limit on 
the parameters. In addition, we assume that both the higgsino and wino mass parameters 
are real and positive. However, as long as one chooses the phases constructively, the results 
would not change since it is possible to include the extra phases of the higgsino and wino 
masses into the redefinition of the phases </?*. 

In Fig. |5l we plot the proton lifetime as functions of the higgsino mass. Here, the 
wino, sfermioro, and color-triplet Higgs masses are set to be M 2 = 3 TeV, Ms = 10 3 TeV, 
and M Hc = 1.0 x 10 16 GeV, respectively. The solid lines correspond to tan/3 = 5, 10,30, 
and 50 from right-top to left-bottom, respectively. Again the shaded region is excluded 
by the current experimental bound, r(p — > K + v) > 4.0 x 10 33 yrs |13j . It is found that 
the lifetime considerably depends on the mass of higgsino as well as the value of tan/3. It 
illustrates that the higgsino contribution is dominant in a wide range of parameter region. 
Indeed, the contribution gets more significant as the higgsino mass is raised up. We also 
show a similar plot in the case of M 2 = 300 GeV and Ms = 100 TeV in Fig. El In this 
case, large tan /3 region is excluded even if the higgsino mass is around 1 TeV, while with 
a rather small value of tan/3 the proton lifetime easily exceeds the experimental limit. 
Anyway, we have found that in the high-scale SUSY scenario, the minimal SUSY SU(5) 



3 To be concrete, we regard Ms as the stop mass, and all of the other sfermion masses are assumed 
to be degenerate with Ms- Generally speaking, stops are lighter than other sfermions, especially those 
of the first and second generations. So, even though one relaxes the degeneration assumption, one ends 
up obtaining a smaller proton decay rate. 
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Figure 5: Lifetime of p — >■ fC + z/ mode as functions of hh- Wino, sfermion and color-triplet 
Higgs masses are set to be M 2 = 3 TeV, M s = 10 3 TeV, and M Hc = 1.0 x 10 16 GeV, 
respectively. Solid lines correspond to tan/3 = 5,10,30, and 50 from right-top to left- 
bottom, respectively. Shaded region is excluded by the current experimental bound, 
r(p -)> K+v) > 4.0 x 10 33 yrs [TJ]. 



GUT is still alive without any conspiracy of suppressing the dimension-five operators. 



5 Conclusions and Discussion 

In this work, we have evaluated the proton decay lifetime via the dimension-five operators 
in the high-scale SUSY scenario. It is found that the higgsino exchanging diagram gives 
rise to the dominant contribution in a wide range of parameter region. After all, we have 
revealed that the proton lifetime may evade the current experimental limit and, thus, the 
minimal SUSY SU(5) GUT is not excluded in the high-scale SUSY scenario. 

In the hh — Ms case, after Mh c being fixed, the proton lifetime depends only on Ms 
and tan/3. In fact, these two parameters are also crucial for the prediction of the Higgs 
boson mass in the high-scale SUSY scenario. Therefore, since now we know that the 
mass of Higgs boson is 126 GeV, it is possible to relate the Ms and tan/3 in the present 
scenario [30TI33] . Further, if the mass spectrum is somehow fixed, we are able to constraint 
Mh c by requiring the gauge coupling unification [SUEOET]. Precise analyses in this 
direction enable us to predict proton decay rate in this scenario, and future experiments 
may examine the prediction. Such kind of model-dependent study is carried out on 
another occasion. 

While the dimension-five proton decay is suppressed by the heavy sfermion masses, the 
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Figure 6: Lifetime of p — > K+v mode as functions of /x#. Wino, sfermion, and color-triplet 
Higgs masses are set to be M 2 = 300 GeV, M s = 10 2 TeV, and M Hc = 1.0 x 10 16 GeV, 
respectively. Solid lines correspond to tan/3 = 3,5,10, and 30 from right-top to left- 
bottom, respectively. Shaded region is excluded by the current experimental bound, 
r{p -* K+v) > 4.0 x 10 33 yrs p3]. 



dimension-six one through the X-boson exchange does not suffer from such a suppression. 
As referred to above, the GUT scale in the high-scale SUSY is slightly lower than the 
ordinary one. Since the dimension-six proton decay lifetime scales as oc M\ with Mx 
the mass of X-boson, it may be significantly enhanced even by a small change in the 
GUT scale. In such a case, the p — > 7r°e + mode may dominate the p — > K+v mode, and 
both of the modes might be searched in future experiments. For instance, the expected 
sensitivities of the Hyper-Kamiokande with ten years exposure [13] are 1.3 x 10 35 and 
2.5 x 10 34 years at 90 % confidence level for the p — > e + 7r° and p — > vK + modesjj 
respectively, which enable us to explore a wide range of parameter region in high-scale 
SUSY models. 

After all, the minimal SUSY SU(5) GUT is still quite promising, and the proton decay 
experiments may reveal the existence of supersymmetry as well as the grand unification. 



Note Added: While this work was being finalized, we noticed the authors in Refs. [55|l63] 
discussed the dimension-five proton decay in a similar context. In Ref. [55], they have just 
shown dimensional analysis to constraint the dimension- five operators, while in Ref. [63] . 
the proton lifetime is examined in supergravity unified models. 



Recent improvements in the analysis of the K + — > 7r + 7r° decay channel may provide a better 
sensitivity for the p — > vK + mode: 3.2 x 10 34 years at 90 % confidence level [62] . 
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Appendix 

A Renormalization Factors 

Here we present the explicit expressions for the renormalization factors defined in Eq. (TT5|) . 
First, we write the renormalization factors A^ and Ar as the products of the long- and 
short- distance renormalization factors: 

Ai,j) _ A Aid) 
A R = A L As , (20) 

where Ai and Al represent the long-distance QCD renormalization factors between the 
electroweak scale (/i = mz with mz the Z-boson mass) and the scale of /x = 2 GeV, 
while Ag and As correspond to the short-distance renormalization effects between the 
electroweak and GUT scales. All of the effects are to be computed at one-loop level. 

A.l Long-range Factors 

First, we discuss the long-distance renormalization factors. The factors consist of two 
effects; one is the running of the quark masses from jjl = 2 GeV to \x = mz, and the 
other is the renormalization effect of the four- Fermi operators in Eq. (115J1 from /j = mz to 
ji — 2 GeV. We neglect the QED corrections since the electromagnetic coupling is much 
smaller than the strong coupling. 

Before analyzing the renormalization effects, we first discuss the input parameters for 
quark masses. In Ref. [64], the values of the light quark masses m q (q = u, d, s) are given 
in the MS scheme at /i ~ 2 GeV, while those for c- and 6-quarks are presented in the MS at 
/i = m c and m&, respectively. Since we define m q in Eq. ( TT5l) in the DR scheme, we convert 
the input parameters into those in the DR scheme. We use the one-loop relation [65] : 

mM = mM{l-^\ , (21) 

where a s = g 2 s /^ with g s the strong coupling constant. For c- and 6-quarks, we evolve 
the masses to \i = 2 GeV by using the RGEs. For top quark, on the other hand, the pole 
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mass is displayed in Ref. [01]. The relation between the pole mass and the DR mass is 
given by 



m t = m t ((j) 



i + ^feiogA + s 

Sir \ m t 



(22) 



Now we consider the QCD renormalization effects. The RGEs for the Wilson coefficients 
of the effective operators in Eq. ( TToT) at one-loop [66] are given aq3 

d Aq 2 

By using the equation, as well as the RGEs for the quark masses, we readily obtain the 
long-range renormalization factors Al and A^. 



At — 


m Ui m di {m z ) ■ C{2 GeV) 


18 18 

fa s (2 GeV)y-s( a s (m b ) y -* 




^■L 


m u Jn di (2 GeV) ■ C{m z ) 


\ a s (m b ) J \a s (m z )J 




A L = 


m 2 m di (m z )-C(2GeV) 
m 2 t m di (2GeV)-C(m z ) " 


fa a (2GeV)y*fa a (m b )\-% 

V a s (m h ) ) \a s {m z )J 


(24) 


Numerically, we 


have 








A L = 0.53, 


A L = 0.34, 


(25) 



at one-loop level. 

A. 2 Short-range Factors 

Next, we evaluate the short- distance renormalization factors. They are composed of three 
factors. First, the effective operators given at the GUT scale receive the renormalization 
effects as they are taken down to the electroweak scale. Second, the Yukawa couplings 
in the color-triplet Higgs exchanging process are determined through the running of the 
couplings from the electroweak scale to the GUT scale. Third, the interaction vertices in 
the one-loop diagrams in Fig. [2] are obtained by evolving the SU(2)^ gauge coupling and 
the Yukawa couplings according to the RGEs from the electroweak scale to the SUSY 
breaking scale, /x = Ms. 

Let us begin with the running of the Yukawa couplings. Initial values for the Yukawa 
coupling constants are given by 

Vuiimz) = J m u .(m z ) , 
\/2mw 

VdX m z) = J m di (m z ) , 
V 2mw 

y*i(mz) = -t= m ei (m z ) . (26) 

V217IW 



Two-loop effects are also calculated in Ref. [57] > though their contribution is found to be small. 
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In the SM, the Yukawa couplings flow according to the following RGEs at one-loop level: 
d 1 



d 
d 

where Y 2 is given as 



167T 2 

1 

16tt 2 

1 
16^2 



yd, 



/ 2 2\ t r 2 2 o2 

2^ ~ 3/dj) + y 2 - ^i - ^2 - 8^3 

/ 2 2\ t 7- 2 2 2 



2^ + ^2 - -<?? - -gl 



^ = E( 3 ^+ 3 4+^), 



(27) 



(28) 



and we neglect the off-diagonal components for simplicity. The equations are also ap- 
plicable when the renormalization scale exceeds the gaugino masses. Above the higgsino 
mass, Eq. (127)1 is valid except that Y 2 is modified to 



^ = £(3^+3^+^ + ^ + ^. 



10* 



(29) 



Here we assume the ordinary supersymmetric relation for the gaugino-Higgs-higgsino cou- 
plings. The couplings may deviate the relation when the SUSY breaking scale is much 
higher than the gaugino and higgsino masses, but it is found that the deviation is usually 
not so significant [231 EM]. 

At the SUSY breaking scale, the Yukawa couplings yj are matched with the super- 
symmetric ones, J/*, as follows: 



VuXMs 



1 



-,yuAM s 



sin/3' 
yd i (M s ) = ZZZ7j y di (M s) 



VeS M S 



cos j3 
1 

cos/3 



VeAMs) 



Above fi = Ms, the RGEs for the Yukawa couplings are given as 



(30) 



d _ 
d _ 
d _ 



Dui 



167T 2 
1 

m^ Vdi 

i 



13 



3 £ Vli + 3 yli + Vl - T79 2 i - 3<7 : 



16 



-9 3 



£(3^ + y%) + 3^ + y\. - —gl - 3g 



16 



-9 3 



IQtt 



\Vei 



E(3^.+<.) + 3^--^ 2 -3^ 



(31) 
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Then, at the GUT scale, the Yukawa couplings f u . and /^ in Eq. §§§ are defined by 

f Ui =y Ui (M Hc ), 

f di =y di (M Hc ). (32) 

For the gauge couplings, the one-loop gauge coupling beta function coefficients are 
given in the SM as 

b a = (41/10, -19/6, -7), (33) 

for U(l)y, SU(2]r, and SU(3)c, respectively. Above the gaugino threshold, the coefficients 
are converted td_| 

b a = (41/10, -11/6, -5), (34) 

and after the higgsinos showing up, they lead to 

b a = (9/2, -7/6, -5) . (35) 

Finally, in the MSSM, they are given as 

b a = (33/5, 1,-3). (36) 

The short-distance renormalization factors for the four-Fermi operators in Eq. ffT5l) 
are presented in Ref. [66]. For the effective operators generated by the wino exchanging 
diagram, the renormalization factor for the Wilson coefficient is 

CU) = ( aM ) ^ ( aM \ "^ ( aM \ '^CM , (37) 

\a 3 ([i )J \a 2 (no)J \ai(no)J 

while for those induced by the higgsino exchange, we have 



' 206-, 



cw=(^t4)^(^t4)" ,2 (^t4) "SO*)- < 38 > 

\a 3 {fXo)J \a 2 {no)J \ai(fio)J 

In the case of \i > Ms, the theory is to be regarded as supersymmetric, and the renor- 
malization factors are obtained as the product of the wave function renormalizations of 
the fields in the effective operators [68J. The wino contribution to the effective operators 
in Eq. ( !T5|) is induced by the effective operators with a form like 

Cij / d Oea^ersetuQ^ 7 Q \ Q\ L™ , (39) 

with i, j = 2, 3, while the higgsino contribution is generated by 

C J d 2 ee a(3 ''U 3a E 3 U ip D h , (40) 

6 When evaluating the renormalization factors, we take the gluino mass equal to wino mass for 
simplicity. 
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Figure 7: Left: A^' j) with 



we set tan/3 



(2,2), (2,3), (3,2), (3, 3) as functions of M s . Here 
3 TeV, and M Hc = 1.0 x 10 16 GeV. Right: A {2 ' 2) as 
functions of Ms- Same parameters as in the left graph are used except for tan/3. Each 
line corresponds to different values of tan/3 (tan/3 = 3, 5, 10, 30, 50). 



L R 

3, M s 



hJ) 



with I = 1,2. Then, the REGs for the Wilson coefficients of the operators are 



d_ 

' dfx 



^'Ei~^' i J 



1 



167T 



[2(f& + vl) + Vi ~ 8g 2 3 - 6g 2 2 - -g 2 ] Q 



V > 



and 






1 



\2{t + t)-^t- l 49l]C. 



(41) 



(42) 



dfi" 16tt 2L wi ' ar ' ~" d 5 

Note that in this case it is important to take the Yukawa interactions into account since 
the effective operators contain the third generation chiral superfields. 

With the RGEs presented above, we finally compute the short-distance renormaliza- 
tion factors as follows: 



.4 



W) 



m 



2 W sin 2/3f u J d . a 2 (M s ) C(m z ) Cy (M s ) 



A, 



An m u .m dj (m z ) a\{m z ) C(M S ) C i:j (M Hc ) 
m A w sin 2 2/3f t y t (M s )y T (M s )f di C(m z ) 



(Aii) 2 al(m z )m 2 t m T m di (m z ) 



C(M^ 



(43) 



H Cl 



With the results obtained above, we compute the renormalization factors. We present 
A^ and A R in Figs. [7] and El respectively, as functions of M s . In both figures, we take 
M s = Hh, M 2 = 3 TeV, and M Hc = 1.0 x 10 16 GeV. In the left graph in Fig. each A { ^ j) 

(2 2) 

is presented with tan p fixed to be tan p = 3, while in the right graph the behavior of A R ' 
is shown for different values of tan /3 (tan (3 = 3, 5, 10, 30, 50). Similarly, each line in Fig. [8] 
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M 2 =3TeV M Hc = 1.0x10 16 GeV 




1(T 10' 10° 
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Figure 8: A R as functions of M s . Here we set M s = Hh, 
1.0 x 10 16 GeV. Each line shows different values of tan/3. 



Mo 



3 TeV, and M Hc 



corresponds to Ar evaluated with various tan/3's. It is found that the renormalization 
factors decrees as the SUSY scale increases, while their dependence on tan /3 is somewhat 
complicated. In addition, the left panel in Fig. [7] illustrates that the effects of the third 
generation Yukawa couplings are significant. 

B Formulae for Proton Decay Rate 

In this section, we display the formulae for the partial decay widths of p — )■ K + v channels 
as well as the hadronic matrix elements which we need to evaluate the decay widths. 
Let us start with the matrix elements. We divide the derivation into two stepsjj First, 
we express them in terms of the low-energy constants a p and /3 P by using the chiral 



7 Calculation of the matrix elements is also conducted by using the direct method j69j . in which 
the three-point correlation functions relevant for proton decay are directly computed on the lattice. 
Preliminary for the recent progress is reported in Ref. [70 . 
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"p ) 



"p ) 



"P ) 



perturbation techniques [7TH75] : 

(K+\e aM {uls{)dl\ P ) = ^-(^-^y LUp 

(K + Mu%di)si\ P ) = -2l-(i + £^.^y LUf 

(k + Mu«A)<p l \p) = ^X^m^) PlUp ' (44) 

where /„. ~ 92.2 MeV [61] is the pion decay constant and the baryon-meson couplings D 
and F are given as D ~ 0.80 and F ~ 0.47, respectively. m p denotes the proton mass, 
while Mb represents the baryon mass parameter in the chiral Lagrangian, which we choose 
as Mb — {m^o + m\o)/2 with m s o and m^o the masses of S° and A , respectively. u p 
is the four-component spinor wave function of proton, and Pl is the projection operator 
defined by Pl = (1 — 7s)/2. The low-energy constants a p and f3 p are defined as 

(0|e a/ g 7 (u^dJX|p) = a p P L u p , 

(0|e^ 7 K4)<b) = P p P L u p , (45) 

with |0) the vacuum state. Second, we determine the constants a p and (3 P . We extract 
them from the results of lattice simulations [74] : 

a p = -0.0112 ± 0.0012 (stat) ± 0.0022 (syst) GeV 3 , 

(3 P = 0.0120 ± 0.0013 (stat) ± 0.0023 (syst) GeV 3 , (46) 

where they are evaluated at /i = 2 GeV. 

With the matrix elements and the effective operators in Eq. (fl5|) . it is straightforward 
to derive the partial decay widths of the p — > K + v^ and p — > i^ + z/ T channels. The result 
is 

r(p -► /f+a) = ^^ICil 2 ( _ m|V ( 47) 

UP ^ A ^ 64 7 r^M^m^sin 2 2/3v mg J ' ^ 

with (i = /i, r) and 

C M = 2/3 p F(M 2 , M 2 ){1 + (D + F)^}fn s V: s J2 muV Uid V UiS e^A% 2) , (48) 



i=2,3 



Win i -~ V^ „ T . .-,„. .(i,3) 



C T = 2/3 p F(M 2 , M|){1 + (D + F)-^}m b V: b ^ m^K^^A 

5 «=2,3 



R 



m 2 m r Kte^i , „,,- f T , Tr / D + 3F m p \ T , Tr 2DmJ 

«p 2 • fl *W Mj ^ J m d K^ s 1 + -/■ + m s K s V^— -/■ . 

77% sin 2/3 I V 3 M B / 3 Mg J 

(49) 
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